Anomalous diffusion and generalized Sparre- Andersen scaling 
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We are discussing long-time, scaling limit for the anomalous diffusion composed of the subordi- 
nated Levy- Wiener process. The limiting anomalous diffusion is in general non-Markov, even in the 
regime, where ensemble averages of a mean-square displacement or quantiles representing the group 
spread of the distribution follow the scaling characteristic for an ordinary stochastic diffusion. To 
discriminate between truly memory-less process and the non-Markov one, we are analyzing deviation 
of the survival probability from the (standard) Sparre- Andersen scaling. 
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I. INTRODUCTION 

The question of the time that it takes for stochastic 
process to reach a specific point or state by the first 
time is central in many applications of stochastic mod- 
eling in physics (Kramers problem [l|), chemistry(reac- 
tion kinetics 0]), biology (neural activity models ]1|) and 
economics (estimation of the ruin time [1]). For a ran- 
dom walk sequence, a nontrivial theorem due to Sparre- 
Andersen 0,0 states that the asymptotic form of the 
probability Sn{x[j) of not crossing the boundary within 
the first n steps after starting the motion at 2:0 > 
(otherwise named the survival probability on the posi- 
tive semi-axis) does not depend on the form of a jump 
length distribution if only it is symmetric, continuous 
and Markovian. For a large number of steps, one in- 
variably has regardless of the exact jump length distri- 
bution type) Sn{xo) oc c(a;o)n~^/^, where the prefactor 
c{xq) depends on the initial position. The result can 
be easily generalized for the continuous-time version of 
the process. For unbiased, continuous Gaussian random 
walk its first passage time density (FPTD) from xq to 
X can be easily calculated explicitly Let us denote 
p{x2, t\x{j, 0) = p{x2 — xq) the probability of motion from 
a position xq to X2 in time t with x denoting the posi- 
tion on the way from xq to X2, i.e. xq < x < X2- By 
taking /(x, r|a;o, 0) = f{x — xq,t) as the probability of 
arriving for the first time at x at time r, the equation for 
p{x2 — XQ,t) reads 

P{X2 - XQ,t) = / p{x2-X,t-T)f{x~XQ,T)dT (1) 



(2) 



with its Laplace transform given by 

P{X2 - Xo,s) = p{x2 - X, s)f{x - Xq, s). 



For unbiased Gaussian random walk we have 
p(x,t|0,0) = p{x,t) = (27rcr2i)-i/2exp(-2f^). With 



the Laplace transform p{x, s) 
{2sa^)~^^'^ exp{— y/^x^sja^) one obtains 



Jo° Pix,t)e '^^dt 



f{x - xq, s) = exp(-v^2(a; - XQ)^s/a^) 



(3) 



which, by inverse Laplace transform, yields the Levy- 
Smirnov distribution 



f{x,t\xQ,0) = f{x 



xo,t) = - 



1 {x-xq)^ 
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(4) 

where xq represents the initial condition. This "in- 
verse Gaussian distribution" decays for long times as 
f{x,t\xo,0) oc and does not have a first moment, 

i.e. the mean first passage time from xq to x diverges. 
On the other hand, since /(x, t|a;o: 0)dt — 1, the par- 
ticle performing the one dimensional Gaussian random 
walk will certainly hit any point x during its motion. 

Assuming the absorbing boundary located at the ori- 
gin, i.e. at a; = 0, formula ^ with x = gives the 
probability density of the first passage time from the pos- 
itive semi-axis for a Gaussian random walk. It should be 
stressed, however that for generally non-Gaussian noises, 
the knowledge of the boundary location may be insuffi- 
cient to specify in full the corresponding conditions for 
absorption or reflection 0, H, H, In particular tra- 
jectories of Levy walks may exhibit discontinuous jumps 
and in a consequence, the location of the boundary itself 
is not hit by the majority of sample trajectories. In order 
to properly take care of possible excursion of the trajecto- 
ries beyond the location of the boundary (at, say x = 0) 
with subsequent re-crossings into the interval {x > 0), 
the whole semi-line {x ^ 0) has to be assumed "absorb- 
ing" . This nonlocal definition of the boundary conditions 
secures proper evaluation of the first passage time distri- 
bution and survival probability 0, llOl ]. see below. 

The very same scenario, see Eq. (U), dictated by the 
Sparre- Andersen theorem holds also true for "paradox- 
ical" diffusion-like processes studied in terms of CTRW 
(continuous time random walks) where kinetics of the 
walker is determined by the distribution of jump lengths 
and distribution of waiting times before a next jump to 
occur pdj . If the process is regular in time but with non- 
trivial jump distribution following the (symmetric) Levy 



2 



law of stability (so called symmetric Levy flight), the 
first passage time density (FPTD) follows the Sparre- 
Andersen universality 0, 0, [13] ■ Notably, however, if 
subordinating the number of steps n to the physical clock 
time t such that the number of steps n per unit of phys- 
ical time follows some distribution with a power-law tail 
p{n,t = 1) oc n~'^^+''^ with < /3 ^ L the deviations 
from the universality can be observed 

To further elucidate the nature of deviation from the 
"standard" Sparre- Andersen scaling in subordinated sce- 
narios, we consider the process X{t) = X{St), for which 
the parental process X{t) is described by a Langevin 
equation 



(5) 



driven by a symmetric a-stable Levy motion La (r) with 
the Fourier transform l^g^^^'^a{T)\^ _ g-T"|fc| _ Here r stands 
for the operational time scale which is changed to the 
physical time scale t by subordination X{St). The sub- 
ordinator St is defined as St = inf {r : U{t) > t} with 
U{t) denoting a strictly increasing i/-stable Levy motion 
(0 < 1/ < 1) and is assumed to be independent from the 
noise term Lair). 

The above setup has been recently proved [H, [13, [lB| 
to give a proper stochastic realization of the random pro- 
cess described otherwise by a fractional diffusion equation 



dpjx, t) 
dt 



d\x\°' 



p{x,t). 



(6) 



Here ^D] ^ denotes the Riemann-Liouville fractional 

with < a ^ 2 



derivative defined as qD] " f{t) = 
TY~^f{T)dT with < J/ < 1 and ^ 



d\x\' 



stands for the Ricsz fractional derivative with the Fourier 
transform T[-^f{x)] = -\k\°'f{k) Occurrence 

of the operator oDl~'^ is due to the heavy-tailed wait- 
ing times between successive jumps and presence of the 



Riesz fractional derivative 



is a consequence of the 



Levy-fiight character of the jumps. 

In this paper, instead of seeking an analytical solu- 
tion to E g . (| 6l) , we switch to a Monte Carlo method 
[3, [I3j [3 I22I [23} which allows generating trajectories 
of the subordinated process X{t) with the parent process 
X{t). Furthermore, we study the potential free case, see 
Eq. ([6]), i.e. we assume V{x)=0. The assumed algorithm 
provides means to investigate the competition between 
subdiffusion (controlled by j/-parameter) and Levy flights 
characterized by a stability index a. 

For Markov processes, the Sparre- Andersen scaling 
[E) [^ presents a universal law which is independent of 
detailed properties of the jump length distribution (if it 
is only continuous and symmetric). In particular, for con- 
tinuous times, the scaling predicts the decay of the 
survival probability, independently of whether the mo- 
ments of the underlying jump process exist or not. For 
example, for a < 2, the moments of the process X{t) (cf. 



Eq. ([5])) exist only for S < a with obvious divergence of 
moments of order S ^ a, i.e. 



(1^1') 



|x|''p(2;, t)dx 



00. 



(7) 



This divergence can be easily demonstrated in the case 
of (pure) Levy flights described by Eq. (O, for which the 
operational time r and physical time t is equivalent, i.e. 
St — t and consequently t = t, see below. In such a case, 
the p{x, t) is a Levy stable density (whose width is grow- 
ing with time) and stays infinite for 5 ^ a. Clearly, 
for finite time t and finite number of representative tra- 
jectories N (otherwise called realizations of the process 
X{t)), variance {X"^) of (symmetric) Levy flights stays 
finite, see [U, Eq. (1.19)] and [l^. In fact, finite number 
of realizations (to be distinguished from the number of 
steps n used in simulation of a single trajectory of time 
duration t = nAt) and finite time introduce an effec- 
tive cutoff to the jump length distribution. In contrast, 
for any fixed time the variance diverges with increasing 
number of simulated trajectories N . Analogously, for any 
fixed N , the variance diverges with increasing time (scal- 
ing like see jH, Eq. (1.19)] and [H). The problem 
of mathematical divergence can be resolved either by in- 
troducing spatiotemporal coupling (typical for so called 
Levy walks [2^. [2^) or by proper truncation of the jump 
length distribution [H, [23, [sO, [U. Several truncation 



methods have been proposed 28l. |29|. [sol [sH to retain the 
finite second moment. In particular, paralleling the sim- 
ulation studies of Mantegna and Stanley [28|] , a smooth 
exponential cutoff has been introduced by Koponen [2^ . 
Instead of truncating tails of a distribution, this approach 
is based upon the exponential tempering of the Levy den- 
sity and preserves the infinite-divisibility [3] of the dis- 
tribution. The classical tempered stable distribution has 
been further generalized by Rosihski_ (for more detailed 
discussion see [H, Chapter 5.7] and j29i. l33l. i34|). 



II. METHODS & RESULTS 

For systems driven by symmetric process the general- 
ized Sparre- Andersen scaling [3, [3, M, [Ij HH] can be used 
to discriminate between Markovian and non-Markovian 
situations. More precisely, according to the Sparre- 
Andersen theorem for a stochastic processes driven by 
any symmetric white noises, the first passage time den- 
sities, f{t) = from the real half line asymptotically 
behave like 



/(i)oct-3/2. 



(8) 



Consequently the survival probability S(t), i.e. the prob- 
ability of finding a particle starting its motion at a;o > 
in the real (positive) half line, scales like 



(9) 
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Therefore, any deviation of the survival probabihty from 
dependence indicates violation of assumptions as- 
suring the proof of the theorem. It can mean either that a 
system is driven by non-symmetric or not "memory-less" 
driving. In consequence, for symmetric drivings, analy- 
sis of data based on (assumed a priori) Sparre- Andersen 
scaling may reveal deviations from Markovianity. 

We study statistical properties of a symmetric free 
Levy motion Eq. ([5]) constrained to the initial position 
a;(0) = 1. To achieve the goals, we use the scheme of 
stochastic subordination [ij, [ij, [11, [2^, i.e. we ob- 
tain the process of primary interest X(t) as a function 
X{t) = X{St) by randomizing the time clock of the pro- 
cess X{t) using a different clock St- The parent process 
X(t) is composed of increments of symmetric a-stable 
motion described in an operational time r and in ev- 
ery jump moment the relation U{St) — t is fulfilled. 
The (inverse-time) subordinator St is (in general) non- 
Markovian hence, as it will be shown, the diffusion pro- 
cess X{St) possesses also some degree of memory. 

The survival probability, see Eq. ([9]), was estimated 
from ensemble of trajectories of the process X{t) start- 
ing at xo {xo > 0). For a < 2, in order to correctly 
account for non-local boundary conditions 0, M, we 
have excluded multiple recrossing events, i.e. every time 
the particle reached any point x beyond the boundary it 
was removed from the system. 

In Figs. the survival probability S{t) = 1 - F{t) is 
depicted for various stability indices a and various sub- 
diffusion parameters v. It is clearly visible that the sur- 
vival probability S{t) behaves like a power-law for all 
considered values of the subdiffusion parameter v and 
stability index a. However, the exponent characterizing 
the power-law dependence is equal to —1/2, as predicted 
by the (standard) Sparre- Andersen theorem, only for the 
Markovian case {v = 1). In more general case the power- 
law is characterized by the exponent b 

S{t) = 1 - F{t) oc (10) 

which differs from 



For a = 2 with any v {{) < v ^1), the first passage time 
distribution is one sided Levy distribution characterized 
by the stability index u/2 |36, i57., j58j, i.e. 



Furthermore, in the general case, the value of the expo- 
nent b does not depend on the stability index a of the 
jump length distribution (35| . 

Figs. [THU confirm that the value of the exponent h de- 
pends on the subdiffusion parameter, v, only. Fig. [2] 
shows results for a = 1.1. Results for others values of sta- 
bility index a are the same as those one for a = 1.1. Fi- 
nally, Fig.[3]presents value of the exponent &, see Eq. (ITUl) . 
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FIG. 1: Sparre Andersen scaling Q for p = 1 (top panel) 
and V — 0.8 (bottom panel) with various a. The process 
was numerically approximated by subordination techniques 
with At = 10~^ and averaged over A'^ = 10^ realizations, 
a = I/V2. 



as the function of the subdiffusion parameter v and sta- 
bility index a. Fig. [3] confirms that exponent b depends 
on the subdiffusion parameter v and the influence of the 
stability index a is negligible. Furthermore, b depends 
linearly on b = -(0.54 ± O.OlW + i0.03 ± 0.03), what 
agrees with earlier findings [H, [H, [sMi , see Fig.[3l Value 
of the exponent 6 is the decreasing function of the subd- 
iffusion parameter v leading to the slowest decay of the 
survival probability for small values of the v parame- 
ter, i.e. when the exponent ly deviates the most from 
its Markovian - "memory-less" value - 1. The devi- 
ation of the exponent b from —1/2 clearly indicates a 
typical slowing down of the subdiffusive process in com- 
parison to its (Markov) regular diffusion analogue. The 
Qf-independence of the survival probability S{t) in this 
case shows that the properties of the decay kinetics are 
determined by the subdiffusive part of the process only. 
This observation is different from the results obtained 
by Sokolov and Metzler for a class of Levy random pro- 
cesses subordinated (via the relation connecting distri- 
bution of number of jumps n in physical time t) to Levy 
flights or to Brownian random walks. In particular, in 
their derivation of subordination, the Authors are us- 
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ing the Markovian Levy flight process X{t) transformed 
to the process X{T{t)) by use of the operational time 
T which, by itself, is called the directing process T{t). 
The density for the process X{T{t)) assumes the form 
P{x,t) = p{x,T)r{T,t)dT -with p{x,t), r{T,t) repre- 
senting densities of a Levy flight process and the density 
of the directing process, respectively. If X{t) is a stable 
process with a stability parameter a and T{t) is a one- 
sided stable process with exponent v, the subordinated 
process X{T{t)) becomes a stable process with the sta- 
bility index va. In contrast, in more general terms of 
the CTRW scenario, after waiving the assumption about 
independent increments of the T{t) process, the asymp- 
totic form of the distribution P{x, t) can be derived by 
use of Tauberian theorems [iGj and is known to be v/a 
self-similar, i.e. = t-''/"'P{xt-''/" ,1) 



in time and superdiffusion in step lengths, the result- 
ing process violates the ergodicity (in the weak sense) 

-0.2 




FIG. 2: Sparre Andersen scaling ^ for a — 1.1 with var- 
ious v. Simulation parameters as in Fig. [l] In Figs, [l] - [2] 
initial position was set to a;(0) = 1. However, due to the 
Sparre-Andersen theorem, results with other values of x(0) 
are perfectly coherent with results for a::(0) — 1 (not shown) 
and lead to the same values of exponent b, see Eq. (|10|l . 



III. SUMMARY & CONCLUSIONS 

We have discussed effects of the subordination scheme 
leading to the fractional diffusion equation Eq. ([6]). By 
use of the Monte Carlo method we have created trajecto- 
ries of the process X{t) — X{St) with St being the inverse 
time a-stable subordinator. Since the St process appears 
as an asymptotic one in the CTRW scheme with heavy- 
tailed waiting time distribution between successive jumps 
and the parental process X(t) is assumed symmetric a- 
stable, the proposed subordination [l^ leads to v/a self- 
similar process whose survival probabilities are governed 
by the stability exponent v. Information gained from the 
analysis of generated trajectories brings around further 
confirmation of non-Markov property of the motion [2^ . 
Moreover, due to the interplay between the subdiffusion 
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FIG. 3: Exponent 6, see Eq. pO|l . characterizing power-law 
behavior of the survival probability S{t) as a function of the 
stability index a (top panel) and subdiffusion parameter u 
(bottom panel). Simulation parameters as in Fig. [T] 



so that the long time average is different from the aver- 
age taken over the ensemble of trajectories [l^, |4l|, [43 |. 
This issue is of special interest in the context of single- 
particle measurements [43j which require analysis of time 
series representative for the motion. In this work we 
demonstrate that subdiffusive and non-Markovian char- 
acter of the motion can be grasped by analyzing survival 
probabilities which deviate from the (standard) Sparre- 
Andersen scaling also in those cases when the ensemble 
averages suggest a Brownian diffusion with v/a = 1/2 
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